Dynamical Invariants for Quantum Control of Four— Level Systems 
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We present a Lie-algebraic classification and detailed construction of the dynamical invariants, 
also known as Lewis-Riesenfeld invariants, of the four-level systems including two-qubit systems 
which are most relevant and sufficiently general for quantum control and computation. These 
invariants not only solve the time-dependent Schrodinger equation of four-level systems exactly but 
also enable the control, and hence quantum computation based on which, of four-level systems fast 
and beyond adiabatic regimes. 



I. INTRODUCTION 

Quantum computation is an emerging discipline in 
which quantum physics is used as a computational re- 
source p], 0]. By making use of states and operations, 
which have no classical counterparts, a quantum com- 
puter is expected to execute some hard tasks for a clas- 
sical computer in a reasonable time. Here operations 
are mostly elements of SU(2"), where n is the number of 
qubits, on which the operations act. It is known that any 
classical logic operation may be realized by a collection 
of the NAND gate. The corresponding "universality the- 
orem" is due to Barenco et al. Q. The theorem claims 
that "any unitary gate can be decomposed into one-qubit 
(i.e., SU(2)) gates and the CNOT gates. In other words, 
the set of one-qubit gates and the CNOT gate are univer- 
sal in gate implementations. In many physical systems, 
implementation of a one-qubit gate is often not hard. It 
may be realized by the Rabi oscillation or the Raman 
transition, for example. In contrast, implementation of 
the CNOT gate can be challenging and its realization is 
sometimes regarded as a milestone for a physical system 
to be a true candidate of a working quantum computer 
Later, it turned out that any SU(4) gate, which en- 
tangles a tensor product state, may serve as an element 
of a universal set of quantum gates with the set of one- 
qubit gates [H, 0]. Important exceptions of two-qubit 
gates that are excluded are the SWAP gate and the "lo- 
cal gates" SU(2) ® SU(2). 

The above observations make the importance of im- 
plementation of an SU(4) gate obvious. Adiabatic two- 
qubit gate implementation is limited in time by the cou- 
pling strength between the two qubits on which the gate 
acts. In NMR, for example, the coupling strength J is on 
the order of 1 ~ 100 Hz, leading to the execution time on 
the order of millisecond to second, while the execution 
time of a one-qubit gate is limited by the strength of the 



RF pulse and it is typically on the order of 10 fj,s for a 
7r-pulsc for a hctcronuclear molecule. It may take longer 
for a homonuclear molecule. By considering this large 
difference between the execution times, two-qubit gates 
often become a bottleneck in shortening the execution 
time of a whole quantum circuit. This is the motivation 
of considering non-adiabatic implementation of nontriv- 
ial SU(4) gates. 

Furthermore, there are many nontrivial quantum algo- 
rithms, such as the Deutsch Algorithm, the Grover al- 
gorithm and the Bernstein- Vazirani algorithm, just to 
name a few, which can be demonstrated with a two- 
qubit system and an SU(4) gate. Execution of these al- 
gorithms with a speed beyond the adiabatic limit shows 
the promising future of a realization of quantum comput- 
ing. 



A. Dynamical Invariants 

An alternative way of constructing solutions of the 
time-dependent Schrodinger equation and obtaining the 
time-evolution operator is by means of the eigenstates 
of an operator / = 1^ Q, which is a dynamical invariant 
of the system with a time-independent expectation value 
(I) • A dynamical invariant (DI) or a Lewis-Riesenfeld 
Invariant (LRI) obeys 

(/) = mwm)) = (m\tfiu\i>(o)) = const., (i) 

where U = U(t; 0) is the time evolution operator. Using 
Heisenberg's equation (in units such that H = 1) we can 
restate this condition as 



= = 
dt 



( 2l 

dt 



(2) 



where the superscript (H) denotes a Heisenberg picture 
operator = U^OU, where O is a Schrodinger picture 
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1 We will use I to denote a dynamical invariant and 1 to denote 
2x2 identity matrix throughout the paper. 
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operator. With the LHS being 0, Eq. (J5J simplifies to 

= ^+i[H,T\, (3) 

which is an equation in the Schrodingcr picture. We will 
refer to this Liouville-von Neumann type equation as the 
DI equation. In terms of the eigenstates of /, \cf> n (t)), 
the general solution of the time-dependent Schrodingcr 
equation can be written as 

|*(t))=^c n e fa »W|^(t)) (4) 

n 

and the time-evolution operator becomes 

C7(t;0) = ^e ia "W|^(t))(^(0)| (5) 

n 

where Lewis-Riesenfeld phase a n (t) is given as Q 

a n {t) = j* (<t> n {s)\ U^- - H{a)j \<j>„(s)) ds. (6) 

We observe that eigenstates of I evolve in the following 
simple form 

e ia »(*) \cf> n (t)> =U(t; 0)|<^ n (0)). (7) 

This passage is transitionless in the eigenbasis of /, and 
is not necessarily adiabatic. The power of DIs goes be- 
yond solving the time-dependent Schrodingcr equation 
however. A major problem in adiabatic quantum con- 
trol is that in many cases the evolution is so slow that 
the system may start decohering. Evolution in the eigen- 
states of /, however, is not restricted by the adiabaticity 
condition and can be made to be fast, a feature which 
caused a recent surge of applications [||-[l6[. DIs have 
been used in the context of quantum field theory [l7| . 
Another attractive feature of DI is that once obtained, 
the time-evolution operator of the system can be con- 
structed from its eigenvectors following Eq. ([5]) (which 
can be otherwise obtained by the direct evaluation of the 
time-ordered integral [lH or through Wei-Norman ex- 
pansion [H-ii]). 

B. Quantum Control Based on DIs 

Quantum control is vital to quantum computation. 
For example, in one of the major quantum computation 
schemes, the AQC [I(||27j], a quantum gate is essentially 
an adiabatic quantum control passage. Nevertheless, in 
general, an adiabatic control of a system is very slow rel- 
ative to the decoherence of the system which may render 
an AQC model over the system impractical. To resolve 
this issue, a quantum c omp utation scheme based on DIs 
has also been proposed [231 . Nonetheless, two limitations 
hinder the quantum control methods and hence the quan- 
tum computation models based on DIs arc not widely 
applicable, due to the three difficulties in obtaining a DI 
of a Hamiltonian. 



In principle, a Hamiltonian of a system is a time- 
dependent operator on the Hilbert space of the system, 
so is a DI associated with the Hamiltonian; they satisfy 
the DI equation. The first two difficulties regard the op- 
erator form of the DIs of a Hamiltonian. A Hamiltonian 
can spawn multiple DIs, as aforementioned; hence, one 
must not only know all possible DIs of a Hamiltonian 
but also select from among them the one pertinent to 
the end of the control, and determine whether the choice 
is optimal. These difficulties worsen dramatically with 
the system's level. Unfortunately, the current literature 
has not been able to cope with these difficulties; instead, 
it usually fumbles a way to certain DIs. Note that in the 
case of two-level, however, a classification of Hamiltoni- 
ans and their DIs does exist [29j]. Consequently, insofar 
the quantum control and computation models based on 
DIs are nigh limited to two-level systems [8l-fl(| . S pec ial 
cases in four-level systems have also been studied |28[ . 

Even if the operator form of a DI is known (it may 
or may not belong to a subalgebra of su(4) as we detail 
in Section [TTJ) it is still hard to obtain the DI in closed- 
form by solving the DI equation, which is in fact a set 
of differential equations, the number of which depends 
on the number of degrees of freedom of the system. Be- 
sides partially causing the first limitation, this third dif- 
ficulty circumscribes how and how much one can control 
a quantum system based on the DIs of the system. In 
fact, because of this limitation, in order to gain the con- 
trol of a quantum system via a DI of the system, a spe- 
cial method — the inversely engineered control (IEC) — 
is needed, which was invented in [l3j for pure states and 
extended to mixed states in Q. 

Although this paper is not about quantum control, so 
as to understand some settings for classifying DIs in the 
main text, let us briefly introduce the gist of quantum 
control, in regard of AQC and IEC. To control a quantum 
system one needs to expose the system in an external 
electromagnetic field that interacts with the system, such 
that the total Hamiltonian of the system and the external 
field evolves temporally in the way that takes the system 
to a designated state at some point of time. A system 
under an AQC remains in an instantaneous eigenstate of 
the total Hamiltonian that must satisfy the adiabaticity 
condition, the violation of which would render the AQC 
of the system impossible. 

In the control based on a DI of a system, however, 
the system does not necessarily follow any instantaneous 
eigenstate of the total Hamiltonian of the system and the 
external field but rather always stays in an instantaneous 
eigenstate of the DI. This idea of control is realized by 
the IEC method. In contrast to the logic that one should 
know a Hamiltonian first to determine its DIs, in IEC, 
one only get hold of the matrix form of a Hamiltonian, 
while leaving its physical parameters unknown, assume 
the matrix form of a DI and postulate some simple func- 
tional forms of its parameters, determine the physical 
parameters of the DI by the boundary conditions that 
comply with the wanted initial and final states of the 
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control, and then use the DI equation to nail down the 
Hamiltonian that bears this DI. In doing so, one avoids 
the difficulty of solving the DI equation directly from a 
known Hamiltonian. Aside of this mathematical advan- 
tage of IEC, a physical advantage is that the control can 
be done regardless of the adiabaticity condition and thus 
serves as a shortcut to adiabatic control. But the costs 
are: 1) the Hamiltonian obtained this way may be un- 
physical or impractical 2) a great deal of control range is 
lost due to the restricted postulates one can make of the 
parameters of the DI. 

In fl3| the method is demonstrated in an su(2) prob- 
lem. The DI is taken to be a function of the form 

/ = — (sin 7 cos f3cr x + sin 7 sin f3a y + cos 70^ ) , (8) 

where 7 and (3 are polynomials in t of order 3 and 4 
respectively, by ansatz. The system is assumed to be 
prepared in a shared cigenstate of / and H initially, and 
due to Eq. ([7]) there will not be any transitions in the 
cigcnbasis of /. This situation is similar to adiabatic 
quantum control where the passage is transitionless in 
Hamiltonian's cigcnbasis, but without the adiabaticity 
condition. The corresponding Hamiltonian is found by 
using the DI equation as 



1 

H =2 



sin j3 



sin ft 



cot 7 cos (3 — $ \ a 



(9) 



The coefficients of the polynomials are fixed by employ- 
ing the population inversion condition at a final time t / , 
and requiring [iJ(0), 7(0)] = [H (t f ) , I (t /)] = which en- 
sures that the initial and final states are eigenstates of 
the Hamiltonian as well, tf can be arbitrarily small, as 
long as it does not violate the time-energy uncertainty 
relation 13, 28], resulting in a passage that can be made 
fast. Unfortunately, a precise implementation of the re- 
quired Hamiltonian is not an easy task. 

One may wonder if in IEC it is still difficult to guess 
the matrix form of a DI of a Hamiltonian. Indeed, it is 
still difficult; however, as IEC has been applied to only 
two-level systems so far 0, the difficulty is negligible, as 
one can always cast the DI in a linear combination of the 
three Pauli matrices, which is simple enough and without 
loss of generality. A four-level system has 15 generators, 
rendering the operator form of DI a lot more difficult to 
guess and resulting set of differential equations are less 
likely to be solved analytically. A companion paper ad- 
dresses the issue by offering a Lie algebraic classification 
of DIs for Hermitian, finite-level systems (30[. To this 
end, we assume in this paper four-level Hamiltonians and 
their DIs live in the Lie algebra su(4) as discussed in [3(| ■ 
Below we list our main results. 



i) We apply the key ideas addressed in [30| to four- 
level systems which arc of particular importance in 
quantum computation. 

ii) We show DIs of four-level systems can be classified 
in terms of the maximal subalgcbras of su(4). 

iii) By means of this classification, we construct DIs for 
a family of four-level system, which are sufficiently 
general for applications. 

iv) Our classification indicates a great reduction of the 
complexity in constructing exact DIs of four-level 
systems. 

v) We list explicit set of differential equations for both 
types, and further discuss the cases which are both 
exactly solvable and physically feasible. 

This paper is organized as follows. In the next section, 
we give a review of the two types of DIs and the adjoint 
DI equation. We move on to the classification of su(4) 
DIs with respect to the subalgebra that spans the set of 
generators in the Hamiltonian in the following section, 
and give detailed analysis for each possible case. Finally 
we summarize our results. Throughout the rest of the pa- 
per, summation over repeated indices is assumed unless 
stated otherwise. 



II. OBTAINING DYNAMICAL INVARIANTS 

We summarize the relevant findings of [30| in a self- 
contained manner from an su(4) point of view in this 
section. A given Hamiltonian can have infinite num- 
ber of DIs. Since I 6 su(4) we can expand it linearly 
in terms of a convenient set of su(4) generators 3 with 
time-dependent coefficients as = g,[ A ^(t)A.;. (For 

brevity, we will drop the superscript denoting the rep- 
resentation for the fundamental representation and set 
g,i = gy^ in what follows). We assume H S su(4), since 
we can always drop the identity element which would oth- 
erwise commute with everything per Eq. ((3|) and result 
in the same equation. Letting Sj be the set that generates 
H and ignoring the coefficients ifijk for a moment (we 
use the convention in which the structure constants obey 



[Xi,Xj] = ifijk^k), due to Eq. J3|) 3 will obey 
3 = [Sj,3]. 



(10) 



2 In [2^1 . where a quantum computation model based on DIs are 
proposed, the control method of four-level systems is essentially 
the same as IEC. 



When we choose to include any generator from the set Sj 
in 3, as a consequence of Cartan decomposition, 3 spans 
the minimal subalgebra © that encloses ft, with f) C 3, 
span(3) = & C su(n). We refer this class as S(uperset)- 
type DI, or for the purposes of this paper equivalently 
S(ubalgebra)-type DI, which has been the focus in the 
literature so far. 

Another class of DIs are generated by a completely 
disjoint set of generators from that of Sj. As we will see 
in Section IIII1 the number of DIs obeying fj n 3 = 0, 
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which wc will refer as D(isjoint)-type DIs, is determined 
by the embeddings of the subalgebras of su(n) , once the 
minimal enclosing subalgcbra of Sj is determined. Clearly 
when & = su(n), there are no D-type DIs. We will see 
that in some cases, such as © = so (5) C su(4), a D-type 
invariant can simplify the problem. 

When wc plug I = giXi into the DI equation, we obtain 

= g + iAg, (11) 

hereafter called the adjoint DI equation, where g = 
(<7i, . . . , gi5) T is a vector made of the expansion coeffi- 
cients of /. At first sight, iA is a 15 x 15 anti-symmetric, 
real matrix whose entries are linear combinations of the 
coefficients in H . By casting Eq. (J3)) into 

;/,A, • <j, II.X, (12) 

and comparing with Eq. (|lip . we see that the A matrix 
can be obtained from the Hamiltonian directly as ad(_ff). 
This means A matrix lives in the adjoint representation 
of su(4), that is A <E ad[su(4)] in general. A obtained 
this way will be a block diagonalizcd 15 x 15 matrix, 
with different blocks representing S-typc and D-types (if 
there is any). Once the adjoint representation is known, 
this serves as a very practical way of obtaining A. 

Due to the DI equation, if / is a valid DI and if a Hermi- 
tian time independent operator A is such that [H, A] = 0, 
then I' = C1/ + C2I + C3A where {c,} are constants is also 
a valid DI. This freedom allows one to work with a "min- 
imal" DI in which the time-independent trace part and 
terms that commute with Sj are dropped. We will ex- 
ploit this freedom in what follows, and work with this 
kind of minimal DIs. We remark that since / obeys the 
same Liouville-von Neumann equation just like the den- 
sity matrix p, it is possible to obtain a density matrix 
from a DI vice and versa, some examples given in [3ll l32l] . 
Traceless part of p is always a valid DI, but the converse 
is not true. A density matrix has trace one and is positive 
semi-definitive whereas a DI is not necessarily, however 
this situation may be remedied using a "non-minimal" 
DI when necessary. 

III. FOUR-LEVEL SYSTEMS 

In this paper, we will be working with a form of four- 
level Hamiltonian which is general enough to cover al- 
most all practical applications 

J.o-i ® o-i + hf Vi ® 1 + hfh ®a l (i = x, y, z)(13) 

with at least one non-zero Ji. Each Ji is a coupling co- 
efficient, and each {h^), is a field acting on the first 



3 By doing that we leave out some systems such as ENDOR in- 
volving anisotropy with electron spin and nuclear spin coupling 



su(3) 0u(l) 


3-4 + 34 + 1() + 80 


so (4) 


(3,1) + (1,3) + (3, 3) 


so (4) ®u(l) 


(1, l)o + (3, l)o + (1, 3)o + (2, 2)2 + (2, 2)_ 2 


so(5) 


5 + 10 



TABLE I. Maximal subalgebras of su(4) in the adjoint repre- 
sentation [34l ]. 



(second) qubit, and is sometimes referred to as a con- 
trol parameter. There could be an extra term 1 (g) 1 in 
the Hamiltonian which we deliberately dropped; it would 
commute with everything else and would not appear in 
the adjoint DI equation. We will choose to drop any gen- 
erator that commutes with all the others for the same 
reason. For purposes of quantum control we prefer to 
work with a resprescntation involving elements of the 
form {<7i £g> 1, 1 (£> Oi,Oi <g> <Jj} which allows us to read 
off the physical content directly. Moreover, this repre- 
sentation is compatible with the spinor representation of 
so(4) and so(5), and will be referred as such from here 
on. Another advantage of this spinor basis is that the 
commutator of two spinor elements is never a linear com- 
bination of more than one generators, which lets us use 
Eq. pU)l smoothly. We will use the shorthand notation 
XY to denote o~ x © o~ y interchangeably throughout the 
paper. 

The A matrix corresponding to the subalgebra (S- 
type) and disjoint (D-type, may or may not exists) cases 
can be read off directly as separate blocks after block- 
diagonalization, in which case g belongs to one of the two 
subspaccs that correspond to two distinct DIs. There 
can be more than one S-type or D-type DIs which will 
cause further block diagonalization in the respective sec- 
tor. The overall block structure and the algebraic struc- 
ture are determined through the embedding of the corre- 
sponding subalgebra of the adjoint representation which 
is the 15 dimensional representation of su(4) in our case. 
Possible group embeddings are listed in Table |U 

In the spinor basis, spanp) can be so(5), so(4) © u(l) 
and so (4) subalgebras of su(4). It can not, however, be 
su(3)ffiu(l) in this representation; one can constraint the 
coefficients of the generators to effectively transform the 
representation into a suitable one but the constraints are 
found to be too tight when we force the "cross-terms" 
such as XY and ZX to be zero in order to adapt Eq. 
(fl~3|) : the number of generators reduces to four, leaving 
us in su(2) © u(l) C su(3) © u(l) (see Appendix E] for 
details). so(4) © u(l) subalgcbra arises from Hamilto- 
nians that can have — but not necessarily present — a 
u(l) generator representing a non-interacting subsystem, 
however because we choose ignore such triviality from the 
beginning we will not have the u(l) generator in the in- 
variant. 

While we can obtain A directly by taking the adjoint 
of H which is practical, we may desire to switch to a dif- 
ferent choice of generator basis afterwards, such as the 
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spinor basis |f| with a particular ordering or a new set 
of their linear combinations. In this case, we may sim- 
ply perform a time-independent similarity transforma- 
tion on the adjoint DI equation: let {Si} be the new 
basis for DI such that / = giXi = g\ The trans- 

formation matrix relating two sets of generators through 
= [SjijXj can be obtained using tv(XiXj) = 25ij as 
[S]ij = tr(Sj Aj)/tr(AjAj). Note that this mixes the orig- 
inal generators rather than transforming them altogether 
in the same way. The g vector will transform the same 
way as A does, since giXi as a whole forms as a single 
"component" of I, from which we deduce that the ad- 
joint DI equation will become 

Q = g^ +iA<P>gV>] 

(14) 

Clearly, for the differential equation 

= ^- t (Sg)+i(SAS- 1 )(Sg) (15) 

to keep its form after the similarity transformation, we 
must have (dS/dt)g = 0, which is always the case when 
the new frame is stationary as in our case. 

We observe that the choice of axes is not of any im- 
portance, as any similarity transformation on the set of 
generators as a whole leaves the subalgcbra spanned by 3 
and the structure constants intact, resulting in the same 
equation as above. Nevertheless, the choice of generator 
basis does make a difference by mixing the components 
of g and determining the block structure of A. 

A. Classification of Hamiltonians and Invariants 

Let us denote an S-type 3 with 3s and a D-type 3 with 
3d- Given a Hamiltonian, the subalgebra spanned by 3s 
can be determined directly by means of Table [ill When 
\Sj\ < 3, due to the simplicity, one can manually check 
whether these generators belong to an su(2). In other 
cases, finding the smallest set in the table such that f) 
fits into gives the answer. However, we can alternatively 
decide on the desired subalgebra as our starting point 
and pick a subset of corresponding generators from the 
table to form a Hamiltonian, which will be our approach 
in this paper to exhaust all possibilities in a concise way. 

Depending on the list of generators appearing in 3s, 
further embeddings are possible through the embedding 
of the subgroups listed in Table [TTJ Note that 50 (4) C 
su(4) and so(4) C So(4) © u(l) C su(4) will result in dif- 
ferent D-types and different block structures for A. The 
block structure is dictated by the maximal embedding. 



su(3) 0u(l) 




so (4) 


{<Ti(Tj, (TjCTj, CTfcl, CTfeCTi, CTfcCTfc, l(Tj}. 


so (4) 0u(l) 


{l(7„cri(7„} U {(Til}, 






{(7il, IcTj, Oj(J k ,(Jj(Ji,<J k Oi,IJ k a k } U {(7, 






{oil, l<Ti, Oj(J k , Ojtjj, o k Oj, a k a k } U {a 


ai}. 


so(5) 


{(Jill, ta*,<Tj<r*,crk<T*}- 



TABLE II. List of maximal (semi-)simple subalgebras & C 
su(4) when fj C {a, 1, la,, a x a x , a y a y , a z a z } with at least 
one <Ji<Ji type generator, in accordance with Eq. (|13|) . There 
is no summation over repeated indices. The indices i,j, k £ 
{x, y, z} are distinct. <r* serves as a wildcard in such a way 
that 0"i<T„ denotes three elements aia x ,<Ji/jy,tJi<j z . The com- 
plementary table is obtained by swapping the role of first and 
second terms in all generators. Alternative tables can be ob- 
tained by similarity transformations. Tensor product symbol 
is omitted for brevity. The table is referred and used in Sec- 
tion ImBi] 

B. Exact Dynamical Invariants 

The adjoint DI equation for the su(2) = so (3) case has 
been studied in the literature extensively under different 
contexts, with many known exactly solvable cases, and 
is known as the Bloch equation with infinite relaxation 
[3oT - [50| . Thus the Bloch equation is a special case of the 
adjoint DI equation. 

In both so(4) = su(2)©su(2) andso(4)©u(l) cases we 
essentially have two decoupled Bloch equations, mean- 
ing that the mathematical problem reduces to solving 
single-qubit problems. Note that the two physical qubits 
are actually coupled in so(4) © u(l) case but not in 
so(4) case. The explicit block diagonalization to bring 
(3,1) + (1,3) into 3 + 3 can be achieved by a time- 
independent similarity transformation. To this purpose, 
we first notice that so (4) is readily split into two dis- 
tinct su(2), {1<7*} U {cr*!} (note that some generators 
may require a correction of sign depending {i,j, k} being 
an odd permutation of {x, y, z} or not). Clearly, such 
a Hamiltonian represents composite systems with two 
non-interacting su(2) subsystems. The other so(4) case 
{liTj, CTfetTj, Ofccr/T U {o-fcl, <Ji<Jj, Cj&j} is equivalent to the 
former case up to a similarity transformation. Moreover, 
the D-type DI for so (4) results in a 9 x 9 A matrix, which 
does not help simplifying the problem. It is because of 
these facts that we will not analyze so(4) subalgebra any 
further 0. In so(4) © u(l) cases however, we respectively 
take the following simple linear combinations of the gen- 
erators to form two separate su(2) algebras 

{(l±<Ti)(7*}/2, 
{<7jt ± lOfc, (TiCTj + CTfcO'i, CfcOj ± CiCi}/2, (16) 
{(T. t l ± lUi^jdj + G),G),.a,(l h ± O-feCTjl/2, 



4 Because there is no standard ordering for the spinor basis and 
we will need to use different orderings later on, we do not list the 
structure constants for them. 



5 Although the physical context is different, [5lll studies an so (4) 
Hamiltonian in a Lie-algebraic framework. 
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in the notation and order of Table [TTJ Note that the sum- 
mands in each new generator are related to each other by 
the u(l) generator of the subalgebra; we have exploited 
the following to obtain these pairs: let {Ti,T 2 ,T 3 } U q 
(where {T lf T 2 , T 3 } C 3 S and q 2 = 1) be an su(2) © u(l), 
then {Ti ± qT 1 ,T 2 ± qT 2 , T 3 ± qT 3 }/2 is a pair of su(2) 
and {T 1 ,T 2 ,T 3 ,qT 1 ,qT 2 ,qT 3 } U {<?} is an so/J) © u(l). 
Different choices such as the one given in (2^, HH can 
be obtained through time-independent similarity trans- 
formations. Furthermore, the sector corresponding to 
the D-type DI further splits into two 4x4 blocks in 
so(4) © u(l), resulting in four smaller, decoupled adjoint 
DI equations. These two sets of 4 generators are conju- 
gate of each other under the respective u(l) generator. 
There are three types of so(4) ffiu(l) (see Table lll|) which 
we will discuss in detail below. 



1. so(4)0u(l) 

1. This first case is also known as single-qubit control 
Ising model, and the corresponding Hamiltonian can be 
written as 



form A (E) = A 



(3,1)+(1,3) W \ U ) W A {2, 2)_ 2 W A (2,2) 2 



H = J XX + h x lX 



hytY 



LIZ. 



(17) 



A quick comparison of the set of generators f) — 
{IX, IF, 1Z, XX} of this Hamiltonian with Tablc|Ij] re- 
veals span(.f)) C so(4) ©u(l) (the first entry in the table, 
due to the distinctive control terms Uer* , along with a^at), 
with the algebraic groupings given as 

{SJ = {IX, 1Y, 1Z, XX, XY, XZ} U {Xt}U 

{ZX, ZY, ZZ, Yl} U {YX, YY, YZ, -Zt] (18) 
= 3 S UQU3 D U3 5 . 

The sets 3d and 3 D representing two D-type invariants 
obey Eq. (fTOjl . and are related to each other by an in- 
finitesimal u(l) charge conjugation 



A)„ where 



^(3,l) + (l,3) 



2i 
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-J 

















(20) 



and (0) is a 1 x 1 matrix containing 0. The corresponding 
linear combinations of so (4) generators listed in Eq. (fl6|) 
are the two independent sets of su(2), {Y>f} = {(1 ± 
X)X, (t±X)Y, (t±X)Z}/2. Transforming to the basis 
= {£+}U{IT} by employing Eq. with [S]^ = 
tr(E^-)/4, 



( 1 1 
10 
10 
10 0-1 
10 

\o o i o 



o \ 


1 



-1 
-lj 



(21) 



will block diagonalize A, 3 i<* + n 3 y yielding 



where 



.4 



-2i 



= g* + iA ± g 



^ — h z h y 

h z ~{h x ±J) 
i — h y h x ± J 



(22) 



(23) 



9 = (Six ±9xx,9tY ± 9xY,9t z ±gxz) T /2. 



[Q, 3 D ] = -2i3 D , [Q,3 D ] = 2i3 D - (19) 

Exploiting the u(l) factor of the subalgebra, it is possible 
to consider an additional time-dependent XI term in the 
Hamiltonian as well, but it will not affect the discussion 
that follows. Although the (3,1) + (1,3) so(4) sector is 
not two readily decoupled su(2) systems, such a situation 
can always be remedied by means of a time-independent 
similarity transformation which we will employ below. 
This is of course not the case for the (2,2)±2 sector. 

In the generator basis the adjoint DI equation is 
= <jr (E) + i^4 (E V S) and A^ has the block diagonal 



2. The second type of Hamiltonian in so(4) © u(l) can 
be written as 

H = JXX + h {1) Yt + h {2) lZ. (24) 

Generators of the S-type, its respective u(l) generator as 
well as the generators of the D-types are grouped as 

{£,,} = {XX, Yl, ZX, ZY, 1Z, XY} U {YZ}U 
{Xt, Zt, YY, YX} U {ZZ, -XZ, IX, —1Y}. [ ' 

Taking the adjoint A^~> = ad(iT), we obtain ApL^j 3 ^ffi 
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hP 






(26) 
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-J 3 


, -hP 
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Similar to the previous case, we note that {XX ± 
ZY, Yt ± tZ, ZX =F XY}/2 generates an su(2) algebra 
which enables us to further block diagonalize -*4.(3,i)+(i,3) 
into 



A d 



2/ 



W±hP 





-J x 

\ v -(4 2) ±4 1) ) J x 



(27) 



with 



9 ± = (gxx ±gzY,g Y t ± 9tz->9zx T .gxy) T /2. (28) 

3. Finally, we consider the Hamiltonian 

H= J X XX + J y YY + h {1) Zt + h (2) tZ. (29) 

which is the last instance of the so(4) © u(l) listed in 
Table M This Hamiltonian can implement a two-qubit 
gate because it can represent Josephson junctions when 
one of the coupling terms Ji vanishes [521; also with the 
addition of the u(l) generator ZZ and the assumption 
Jx = Jy = J zi the Hamiltonian can implement D3 model 
[53l [HiJ and quantum dots [H| . 



The generators are grouped as 

{£,:} = {XX, YX, Zl, YY, XY, 1Z} U {ZZ}U 
{YZ, tY, XZ, IX} U {Xt, ZX, -Yt, -ZY}, 



(30) 



and is A 



.4 



(3,1) + (1,3) 



= 2i 
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(31) 
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{XX T FF,YX ± XY,Z1 ± tZ}/2 is an su(2) and 
j4(3,i)+(i,3) decomposes into 

-(/ l (2) T / l (i)) 
A ± = 2i\ ftP) =p ftP) 



T 




( ± J,) 







9 = (gxxTgYY,gYx±gxY,g z t±9iLz)/ 2 - ( 32 ) 

2. 50(5) 

We now turn to the most general so (5) case of the 
Hamiltonian Eq. (|13| . which also is the last case in Table 

m 



H =J X XX + J y YY+ 

h (2) tX + hf>tY + h (2) tZ + h^Zt 



(33) 



This time, the generators are grouped into two sets per 
Table |U 

{£} = {tX, tY, tZ, XX, XY, XZ, YX, YY,YZ,Zt} U 
{ZX, ZY, ZZ, XI, Yt} (34) 

A^ matrix is a direct sum of two blocks A^ © A^ , 
where 
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This is an example case where a D-type DI greatly 
simplifies the problem. We further observe that there 
are two blocks in A 5 which are "coupled" to each other 
through J x and J y . The first block is a Bloch equation 
similar to Eq. (j2"3)) , whereas the second block is solved 
by 



9x1 + alt = C > 9x1 + 9xi = Ch£ 



(36) 



This feature can be exploited in construction of pertur- 
bativc solutions in the weak-coupling limit. 



C. Solutions of the Bloch Equation 

We have encountered A € so(3) for 3s G so(4) © u(l) 
in Eqs. (J23J) , ([2TJ) and ((32j) . Such a matrix lives in the 
adjoint representation of su(2), as a result we can find a 
two-level Hamiltonian satisfying = ad(_ff ± ), meaning 
we have effectively reduced four-level problems with H £ 
so(4) and so(4) © u(l) into a pair of independent two- 
level problems, associated with the Hamiltonian H = 
(h x ± J)<j x + hyOy + h z a z for Eq. ([23]). (see Table III). 

For the special case of Eq. ([23| with h x = const., 
h y = B cos tut, h z = Bsmtut where B and u> are 
time-independent, the Hamiltonian Eq. ()1T() becomes 
the NMR Hamiltonian in rotating-frame, and the prob- 
lem has a simple exact solution = (±J + h x — 
oj/2,Bcosut,BsmLut) T such that we obtain 



Is =(±J + [hx - w/2])([l ± X]X)+ 

B cos ut([l ± X]Y) + B sin wt([l ± X] Z) 



(37) 



We note that the Hamiltonian and the invariant are both 
periodic with T = 2it/uj and thus can be used to obtain 
cyclic non-adiabatic geometric phases [56l [57| ■ The un- 



normalized eigenvectors and eigenvalues of a simpler DI 
I s =l£ +Ig = 2JXX + (2h x - w)ll-f 



2BcoscottY + 2BsmuttZ 



(38) 



follow as 



= *(gt + *9t) 
A| =±g + ; 

1^) = *(.9 2 ~ +i9i) 
Xt = ±g~ 
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+ (gt±9 + ) 
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(5 3 ± 9' 



( \ 

-1 


V 1 / 



(39) 



where <7 ± = \g |. Using these eigenvectors in Eq. ([5]), 
we can obtain the analytic time-evolution operator of the 
system. It would be straightforward to design quantum 
control of such systems accordingly. 

Finally, the type of Bloch equations that appeared in 
Eq. (f27j) and Eq. ([32]) has several exact solutions [43| 
corresponding to experimentally realizable fields. When 
the time-dependent term is periodic and is much greater 
that the time-independent part (weak-coupling limit), a 
perturbative solution can be constructed j44|- 



IV. DISCUSSION AND OUTLOOK 

In this work we offered a Lie-algebraic classification 
and construction of DIs of four-level systems. Such sys- 
tems cover two-qubit systems which are relevant and es- 
sential for quantum computation. Applications of DIs 
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TABLE IV. A matrices, obeying tr(AiAj) = 28%j. Ai_3 is an su(2) and Ai^g is an su(3). Ag commutes with A1-3, A15 commutes 
with Ai_8, giving rise to su(2) © u(l) and su(3) © u(l) subalgebras. 



&d(a x ) = 2i 



ad(<r z 



2i 
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0-1 

Voi 

/ -1 

1 
Vo 



eA(py) — 2i 



I 1 


\ -1 



TABLE III. Adjoint representation of su(2), [ad((Tj)]jfe 
-i(fi)jk where [<Ji,<Tj] = ifijk^k- 



go beyond the exact solutions of the time-dependent 
Schrodinger equation, such as fast control of four-level 
systems beyond adiabatic regimes and non-abelian geo- 
metric phases. We have shown that the span(3) = so (4) 
and span(3) = so(4) © u(l) cases are reduced to the 
solutions of two-level systems, with analytic solutions 
corresponding to practically realizable four-level param- 
eters. One such solution was explicitly given. The 
span(3) = so (5) was noted for its D-type DI, whose solu- 
tion can be reduced to solving two weakly coupled sets of 
differential equations. These DIs can be applied to gain 



various control of four-level systems. In an ongoing work 
to be reported elsewhere, we try to devise a CNOT gate 
as a control passage based on the DIs constructed in this 
paper. 
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Appendix A: No su(3) S-Type DIs 



J 



The fundamental representation of su(4) is 15 Hermi- 
tian set of {A;} matrices (Tables IWl and N\. The first 8 
of these, along with A15 which commutes with the rest 
form the su(3) © u(l) subalgebra. They are related to 
spinor generators as 

2Ai = (1 + Z)X, 2X2 = (1 + Z)Y, 

2A 3 = (1 + Z)Z, 2A 4 = Y(1 + Z), 

2A 5 =Y(1 + Z), 2A 6 = XX + YY, 

2A 7 = YX — XY, 2V3A 8 = 2Z1 - 1Z + ZZ, 
V&\i5 = Zt + tZ - ZZ. (Al) 



2H written in terms of Xj is 
h^X 1 + h^X 2 + {h^ + J Z )X 3 + 

h { x 1) X 4 + h { y 1 '>X 5 + {J x + J y )Xe + 



2fej 1} - + J z 



A 8 + 



(J x - J V )X 9 + h<p\ u + hW\ 12 + A 13 + ^ 2) A 14 + 
2 



(hW+h?-J z )\ u 



(A2) 



We see that in order to get rid of Ag_i4 terms, we must 
have hx = hjp = h x 2 ^ = hy 2 ^ = and J x = J y . This, 
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TABLE V. Structure constants for A matrices. 



however, reduces the Hamiltonian to 

(/4 2) + J 2 )A 3 + 2JA 6 + lhz -| +Jz X 8 + 
2=(hP+hW- J Z )X 15 , (A3) 

where we defined J = = J y . The coefficients are 
linearly independent but the set of generators belong to 
su(2)©u(l) C su(3)®u(l), from which we conclude that 
one cannot emulate and gain the full SU(3) control of a 
qutrit using two qubits with a Hamiltonian of the form 
Eq. (P3). 
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